Resonance theory of the crossover from Bardeen-Cooper-Schrieffer superfluidity to 
Bose-Einstein condensation in a dilute Fermi gas 
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We present a description of the behavior of a superfiuid gas of fermions in the presence of a 
Feshbach resonance over the complete range of magnetic field detunings. Starting from a resonance 
Hamiltonian, we exploit a functional method to describe the continuous behavior from Bardeen- 
Cooper-Schrieffer to Bose-Einstein condensation type superfluidity. Our results show an ability for 
a resonance system to exhibit a high critical temperature comparable to the Fermi temperature. 
The results are derived in a manner that is shown to be consistent with the underlying microscopic 
scattering physics. 
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I. INTRODUCTION 

The ability to cool a gas of fermionic atoms well into 
the regime of quantum degeneracy hints at the exciting 
possibility of allowing one to study the mechanisms of 
superfluidity in an entirely new context JlJ. Since these 
systems remain extremely dilute and are not complicated 
by long-range Coulomb interactions or lattice effects, cold 
degenerate gases seem ideal for the study of the funda- 
mental physics behind the exotic behavior of superfluid- 
ity. Unfortunately, temperatures of around 0.2Tp are the 
current state of the art in cooling [Q, || |, §, |, . To ob- 
tain such high critical temperatures a strong coupling 
mechanism is required forcing the theoretical descrip- 
tion to extend beyond the standard Bardeen-Cooper- 
Schrieffer (BCS) approach. 

Several theoretical papers have studied the effects of 
increasing the two-particle interactions [@, ||, [)), charac- 
terizing the coupling processes by large negative scatter- 
ing lengths. We have focused in detail upon a related yet 
distinct approach [jlo|| , which is to significantly increase 
the interatomic couplings by making use of a Feshbach 
resonance. The difference is that, in the neighborhood 
of the resonance, the interactions can no longer be ade- 
quately described by a scattering length, since the scat- 
tering length diverges as one approaches the resonance. 
This is an artifact of the approximations made in formu- 
lating the theory since the full energy-dependent scat- 
tering T-matrix, which is the true descriptor of the two- 
particle interactions, does not diverge at finite scattering 
energy. This leads us to explicitly incorporate the physics 
of the resonance into our microscopic description of the 
interatomic couplings. A much more detailed discussion 
of this can be found in Ref. |]l0| . 

So far, we have not discussed in detail the role of fluc- 
tuations, which can have a significant effect on the criti- 
cal behavior jllj . How we incorporate these fluctuations 
proves crucial in describing the physics correctly within 
the crossover regime where we find a significant popu- 
lation of tightly bound composite particles. The aim, 
therefore, of this paper will be to account for fluctuations 
in such a way as to properly describe the behavior of a 



superfiuid Fermi gas at all detunings from the resonance. 
Recently, a complimentary treatment was independently 
developed by Ohashi and Griffin |l3)] . The slight quan- 
titative differences between their results for the critical 
superfiuid transition temperature and the values we will 
present here appear to arise primarily from the use of 
quite different two-body scattering parameters and a dis- 
tinct renormalization procedure. 

The problem of describing a superfiuid Fermi gas at 
all coupling strengths has been extensively studied in re- 
cent years, motivated by a desire to explain the proper- 
ties of "exotic" high-T c superconductors, whose behav- 
ior seems to lie in a region somewhere between BCS 
and Bose-Einstein condensation (BEC) superconductiv- 
ity. An early description of the crossover from BCS to 
BEC superconductivity was put forth by Nozieres and 
Schmitt-Rink (NSR) after the pioneering work of 

Eagles Jl5| and Leggett [fBL and later expanded upon 
by various authors |l7|, |18|, |1S(|. A functional analysis of 
the crossover behavior, which is the method that we will 
employ, was equated to the NSR method by Randeria et 
al. [[L9| . We will adapt this method to a resonant system 
as necessary to describe the relevant physics of superflu- 
idity in dilute atomic gases. It should be stressed that our 
method contains the multichannel interatomic couplings 
intrinsic to the Feshbach resonance. This was not consid- 
ered in previous calculations in the context of condensed 
matter systems. 



II. RESONANT ACTION 

We consider the Feshbach resonance [^O) for s-wave 
scattering of atoms in the lowest two hyperfine states 
of a fermionic alkali atom, denoted symbolically by 
a £ {t, !}• For a homogeneous system we have the fol- 
lowing generalized Hamiltonian: 

H(t) = E / ri(x)(H«-v)Mx)d 3 x (1) 
+ J2 J ^ mj (x)(H mj -2n + v J )y) mj (x)d 3 x 
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U(x - x')^\(x)ijl(x')^ l (x')^(x)d 3 xd 3 x' 



J 



-^^x^ix') +H.c. d 3 xd 3 x', 



where the operators create (annihilate) fermions 

at x — (x,t), and ipjn- Ww) create (annihilate) compos- 
ite bosons. The free dispersion Hamiltonian for fermions 
(bosons) is H a (H mj ) and Vj is the detuning of the j th 
molecular state from the collision continuum. The col- 
lisional interactions are described by both background 
fermion scattering (U) and an interconversion between 
composite bosons and fermion pairs (gj). 

Functional methods prove to be especially convenient 
in describing the thermodynamics of the resonant system. 
For a finite temperature field theory, the connection with 
statistical mechanics is made by Wick rotating the time 
coordinate t — > — it so that one works in terms of the spa- 
tial coordinate x and temperature r 21. In this space, 
we define the action in the usual way 



dr I d xipJ(x,T)d T ij)i(x,T) 







H(r)dT, 



(2) 



where the sum in I runs over both the Fermi and the 
Bose degrees of freedom. In this functional formulation 
we treat the fermion fields Vv as Grassmann variables 
p2| and the composite Bose fields ipm, as classical fields. 

Let us consider a system comprised of fermions at 
some finite temperature r inside a box of volume V 
(for convenience, let us work in the set of units where 
h = kf, = 1). By imposing periodic boundary conditions 
upon the fields ip a and ip m , we form the following Fourier 
series expansions 

W(*,r) = (/^-V^e^+P-^p), 

k.uj 

^ mj {x,r) = (pvr 1/2 Y, el{VT+q ' X)b ^> ( 3 ) 

q . v 

with even thermal (Matsubara) frequencies for the 
bosons (v = 2im//3, where n is an integer) and odd fre- 
quencies for the fermions (u = 2i:{n + l)/(3), to pre- 
serve the particle statistics. Here a cr (p) annihilates a 
fermion at p = (k,u>) and bj(q) annihilates a molecule 
at q= (q,uj). 

By making use of the above transformation, Eq. (^), 
we may write out the action for the resonant system in 
terms of the Fourier coefficients a CT (p) and bj(q). In or- 
der to help clarify the following calculation, we split the 
resulting resonant action into two parts, the first being 
the usual BCS action: 



S 



BCS 



(4) 



ijy £ t/aT(P 1 ) a *(P 2 ) a l(P3) a T(P 4 )- 

Pl+P2=P3+P4 



The remaining part of the action we will label the molec- 
ular action 



~M 



E 1 

1,3 



^ - — - ^ + 2{j)b* j (q)b j (q L ) 



(5) 



mff E 9j ( & j*(q)ai(Pi) a T(P2) + a UP2)a* l (p 1 )b*{q)) 



q=pi+P2 ,j 



In deriving Eqs. (0) and (||) we have inserted contact 
potentials for the couplings U(x — x') — > US(x — x') and 
gj(x — x') — > gjS(x — x'). The full partition function 
for our resonant system, under the model Hamiltonian of 
Eq. (0) can now be written as 



z = j (nixDa^ fn m i D6 ij ° 



Sbcs+Sm 



> (6) 



with the functional integral, Dc = Oi^ *' ranging over 
all Fermi and Bose fields. 



III. SADDLE-POINT APPROXIMATION 

From the form of the action in Eq. (||), it should be 
apparent that all of the resonant contributions are con- 
tained within the molecular action. In practice this gives 
rise to the integral of a displaced Gaussian that can be 
easily evaluated. After integrating out the molecular de- 
grees of freedom, we are left with the partition function: 



z = \ Ii Zs M/ 4m 




(7) 

Here 2^(5? /4m + Vj — 2fi) is a Bose partition function 
describing the formation of bound molecules and SbcS' 
is the BCS action with a potential that is now dependent 
on both thermal frequencies and momentum. The inter- 
action potential in the BCS action is, therefore, modified 
in the presence of a Feshbach resonance in the following 
way: 



U^U- 



q^/^m + Vj - 2fi 



(8) 



With the above partition function, Eq. (]?]), we may go 
on to calculate all thermodynamic properties of interest. 
Here, we are primarily interested in calculating the criti- 
cal temperature of the supcrfluid phase transition. This 
can be done by solving for the gap and number equation, 
and then self-consistently solving these two equations for 
both the chemical potential and the critical temperature. 
The procedure is straightforward since the full resonant 
calculation has been reduced to the usual BCS calcula- 
tion, only with a more complicated potential. Following 
Popov's derivation p2[ , introducing the complex auxil- 
iary Bose field c(q) and expanding about the neighbor- 
hood of its zero value (which is equivalent to saying that 
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we expand about the zero of the gap near T c ), we derive 
the gap equation at the critical point 



2/x 



E 



tanh(/3(fc 2 /2TO-/i)/2)' 
2(fc 2 /2m- y) 



X(q) 



(9) 

The second self-consistent equation, the number equa- 
tion, is found in the saddle-point approximation by ex- 
panding the action to lowest order, i.e. c(q) = c*(q) = 0, 
and using the thermodynamic identity N — —dlnZ/d/j, 
giving 



N 



i 



0,k 



,/3(fc 2 /im+Uj-2^) _ 



1 



e /3(fc 2 /2m- M ) _|_ J" 



(10) 



Thus, our number equation counts all free fermions, Nf, 
plus an additional boson population Nf,. Equations (g) 
and ( [To| ) provide for us a set of equations for determin- 
ing T c and /it at the critical point. In the usual BCS 
theory, this level of approximation proves reasonable for 
calculating T c in the BCS limit (small negative scatter- 
ing length) but diverges as the scattering length grows, 
and is wholly inapplicable for positive scattering lengths. 
The reason for this is that the primary mechanism for 
the phase transition within the weak coupling BCS limit 
is the formation and disassociation of Cooper pairs. As 
the coupling increases the particles tend to pair up at 
higher and higher temperatures which means that the 
critical transition is no longer signaled by the forma- 
tion of Cooper pairs, but rather by a coherence across 
the sample caused by condensation of preformed Cooper 
pairs. Since we are interested in describing the resonant 
system at all detunings, the equations that we have de- 
rived so far are insufficient because they do not account 
for this process. We should, therefore, next focus on 
how to more accurately incorporate atom pairing into 
our model. 



IV. BEYOND THE SADDLE-POINT 
APPROXIMATION 

To account for fluctuations in the fermion field, we 
follow the method of Nozieres and Schmitt-Rink Ej] in 
its functional form as put forth by Randeria et a~19[. 
This procedure will introduce a next order correction to 
the saddle point calculation of the previous section. By 
expanding the action to second order in c(q), and then 
calculating the number equation in the same way as was 
done when deriving Eq. (10), we introduce an additional 
population into the equation. The action becomes 



S(c(q),c*(q)) 



BCS' 



s 



BCS 



(0,0)- 



q 



|c(q)| 2 x(q), (11) 



where we have defined the auxiliary function x(l) as 



U ■ 



■S= + Vj-2fi- iuj 



E 



1 - /(eg+fc) - f(ei-k) 
+ £|_ fc - iuj 

(12) 



Here f(tk) is the Fermi distribution function and 
£fc = k 2 /2m — [i. The resulting modified number equa- 
tion is 



n = 2y: 



i 



- + 2 V- 

e /3(fc 2 /4m+^-2^) _ I jL^i g/3(fe 2 /2m-^) _|_ \ 

j,k k 



1 



£E|;k8[i-x(«,fc>)]. 

q,uj r 



(13) 



This inclusion of the first order fluctuations introduces a 
population of atom pairs, N p , that behave like bosons. 
We are now able to solve for the fluctuation corrected 
critical temperature from a self-consistent solution of 
Eqs. (|) and @. 

Due to the contact form of the couplings that we have 
chosen, however, we are immediately plagued with prob- 
lems of divergences in our equations. This can be reme- 
died by a proper renormalization, which means replacing 
the 'bare' couplings and detunings by the correct renor- 
malized forms depending on the actual physical parame- 
ters as well as a momentum cut-off K cut . This same pro- 
cedure is needed to renormalize the usual, nonresonant 
BCS theory ]2]J] and is done by relating the bare poten- 
tial to the two-particle scattering matrix, T, through the 
Lippman-Schwinger equation. This leads to a renormal- 
ization of the coupling as U = TU, where the bar nota- 
tion represents a bare parameter. We make the definition 
a = mK cut /(2ir 2 ) and introduce the dimensionless factor 
r = (1 — aU)^ 1 . This renormalization incorporates the 
true microscopic physics into the problem, removing the 
unphysical divergence. 

For the resonant case, a more sophisticated approach 
is required, which we have extensively studied previ- 
ously jl(j. In the case of a single resonance (i.e., the 
j sum has only a single term) the renormalization is per- 
formed by the following relations: U — TU, g — Tg, and 
v — 9 + agg. The chemical potential of each atomic state 
is modified to /i = p,— < T(k)nk > to include the proper 
mean-field shifts induced by all two-body scattering pro- 
cesses, where is the Fermi distribution, T(k) is the 
two-body scattering matrix for the resonant system [jl0| , 
and <> denotes an averaging. This shift, however, is 
sufficiently small to neglect; inclusion has demonstrated 
corrections of the order of 1% or less. By replacing the 
bare values in Eqs. and ( |l3"| ) with the renormalized 
values, all of the results to be discussed have been shown 
to be independent of the introduced momentum cutoff 

Kcut- 

The renormalization of the resonance theory forces us 
to take a closer look at the bound state physics of the sys- 
tem. In Fig. (Q) we show the bound state energies for a 
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FIG. 1: Binding energies for 40 K resonance at positive back- 
ground scattering length ab g = 176ao. A single resonance 
with a positive background scattering length produces an ef- 
fective scattering length a c ff = ab g (l — k/v), as seen in the 
figure inset where we plot the effective scattering length vs. 
detuning (the dotted line is at 176oo). A positive a^ g , which 
is larger than the range of the potential, implies that another 
bound state is not far below threshold (dash-dotted line). In 
combination with the Feshbach state (dashed detuning line) 
this results in an avoided crossing and the molecular state of 
interest asymptotes quickly to the dash-dotted line. 
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FIG. 2: Same as Fig. [jj but now for an artifical situation 
with abg < 0. A resonance system with a negative background 
scattering length has only one bound state relatively close 
to threshold, which is shifted positive of the detuning. The 
next bound state in the potential is too far away to be of an 
significant influence. The line styles are the same as for Fig 



we solve the gap equation for the critical temperature. 
The result is the usual exponential dependence on the 
effective scattering length 



single resonance system with a positive background scat- 
tering length. The figure results from a coupled square 
well calculation of the bound state energies p0| , and 
shows the avoided crossing of two molecular states. The 
upper state behaves to a fairly good approximation as 
Eb = (miijj)" 1 , which is the molecular binding energy 
regularly associated with a contact interaction p3| . The 
lower state, however, is offset from the detuning by an 
energy ~ k and goes linear with the detuning. We find a 
similar behavior as in the lower state in the first term of 
Eq. Taking the cutoff to infinity, which is justified 

since this term does not diverge, the renormalized detun- 
ing approaches v — > v — g 2 /U . This produces a constant 
shift of g 2 /U — k between the detuning and the molec- 
ular binding energy. Keeping this term in the number 
equation would incorrectly cause a transfer of the entire 
population into the wrong molecular state. In order to 
avoid this unwanted behavior we set this term to zero, 
i.e. Nb = 0. In the case of a negative background scat- 
tering length, we would not have encountered this prob- 
lem, and only one molecular state would have appeared 
(see Fig. |]). We will show in the next section that the 
pairing term fully accounts for the correct population of 
molecules in this system with a b g > 0. 

Before we present the full crossover solution for the 
case of 40 K, let us look at the analytical solutions to 
Eqs. (|) and (H) in the strong (BEC) and weak (BCS) 
coupling regimes. We will first turn our attention to the 
weak coupling (BCS) regime. In this limit we would ex- 
pect only free fermions to contribute to the population, 
so from equation ([l3]) we find that the chemical potential 
is at the Fermi surface (li = E-p). With this information, 



T C /T F 



=7-2 



exp( 



2fc f |a c ff 



(14) 



where 7 ~ 0.5772 is the Euler-Mascheroni constant, fcf 
is the Fermi wave number, and a c ff < is the effective 
scattering length produced by the Feshbach resonance 
a c ff = db g (l - k/v). 

The other limit we may consider is the strong coupling 
(BEC) limit. When the argument of the tanh function in 
the gap equation (^) becomes sufficiently negative, it is a 
good approximation to use its asymptotic value of unity. 
What this means physically is that the fermion statis- 
tics are unimportant in determining the value of the gap. 
This allows us to solve the gap equation for the chemi- 
cal potential as a function of detuning. In the limit of 
large negative detuning we find that fx — > —Eh/2, where 
Eb ~ l/ma 2 s . Within this limit the entire population 
has been converted to molecules and we can solve the 
number equation to get the BEC condensation tempera- 
ture of T C /T F ~ 0.218. 



V. NUMERICAL RESULTS 

To study the transition between the BEC and BCS 
regimes, we numerically solve Eqs. (0) and (13) for 40 K. 
The single resonance curve is produced using a back- 
ground scattering length of 176ao and K = 7.68G at a 
density of 10 14 cm~ 3 . 

Figure [| shows the critical temperature as a function 
of magnetic field detuning. The crossover calculation 
clearly merges with the BEC result for large-positive 
detunings and smoothly connects between positive and 
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FIG. 3: Critical temperature T/Tf as a function of detuning 
v in gauss. The dashed line corresponds to the usual BCS 
solution, which limits to the full crossover theory at large 
positive detuning. At negative detuning, T c drops to the BEC 
condensation temperature of T c /Tf ~ 0.218. 
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FIG. 4: Chemical potential as a function of detuning v in 
gauss. For large negative detuning 2/j, approaches the bound 
state energy of the molecular state. At increasing positive 
detuning, the chemical potential slowly approaches the Fermi 
energy. 
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FIG. 5: The fraction of the total population as a function of 
detuning v in gauss. The dot-dashed line corresponds to the 
pairing fraction N p and the solid to the free fermion fraction 
N f . 



negative detunings, limiting to the Bose condensation 
temperature of T C /T F ~ 0.218 for large negative detun- 
ing. This approach gives a maximum near zero detuning 
(T c /Tp ~ 0.26) , but the maximum critical temperature 
we find is less than the predictions of the HFB approach 
in our earlier papers (T c /Tp ~ 0.5) [Q. We believe this 
is due to the inclusion of fluctuations in the beyond sad- 
dle point approximation which act to reduce the gap at 
zero temperature and therefore the critical temperature 
for the formation of a superfluid. 

Figure || shows the chemical potential as a function of 
detuning, beginning at the Fermi energy for positive de- 
tuning and approaching half the bound state energy at 
large negative detuning — > —Eb/2. Figure (||) shows 
the change in population as a function of detuning. For 
large positive detuning, the system is composed solely of 
free fermions. As the detuning is decreased (i.e. from 
positive to negative) the contribution of the fermions be- 
gins to decrease until all the population is transferred into 
the atom pairs at v ~ — 0.5G. The chemical potential is 
then equal to —Eb/2 and we may identify the atom pairs 
from that point as the molecules. The superfluid behav- 
ior then comes from the condensation of these molecules, 
which are no longer disassociating into free fermions. 



VI. CONCLUSION 

We have presented a crossover model to describe the 
behavior of a gas of fcrmionic atoms for all detun- 
ings from a Feshbach resonance. The model is able to 
smoothly connect between the BCS and BEC regimes 
and accounts for the microscopic two-body physics 
throughout. We find a smooth behavior of the transi- 
tion temperature in the entire crossover regime, with a 
maximum near zero detuning of T c /Tf ~ 0.26, and agree- 
ing with the appropriate BEC and BCS behaviors in the 
respective large detuning limits. The maximum is be- 
low the value predicted by the Hartree-Fock-Bogoliubov 
theory for the uniform gas derived in earlier work. This 
result is a direct indication of the important role of pre- 
formed atom pairs which are neglected above T c in the 
Hartree-Fock-Bogoliubov theory. Such atomic pairing is 
represented as fluctuations in the fermion pairing field 
and modify the elementary excitation spectrum even in 
the normal phase. We emphasize, however, that in this 
paper we have only accounted for the pairing physics by 
including the second order fluctuations in order to be able 
to account for the correct molecular binding energy as de- 
rived from the two-body scattering physics. While the in- 
clusion of this order of fluctuations is the main ingredient 
necessary to encapsulate essential aspects of the behavior 
of the system in the crossover regime, it would be inter- 
esting to extend the approach to consider the effect of 
higher order interactions which have not been accounted 
for. For example, we have only included interactions be- 
tween free fermion atoms, neglecting all other contribu- 
tions such as the interactions between pairs. This is most 
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clearly seen in the BEC limit, where our solution adopts 
the thermodynamic behavior of the ideal Bose gas, rather 
than the dilute interacting Bose gas. A more sophisti- 
cated treatment could extend our approach to consider 
all these factors. Nevertheless, the results we have pre- 
sented here illustrate the realistic potential in this realiz- 
able system for increasing the superfluid transition tem- 
perature, with the aid of the Feshbach resonance, to a 
significant fraction of the Fermi temperature. This is an 
important and timely aspect from a practical perspective, 
because the maximum value we predict is in the region of 
the temperature range which is currently experimentally 



accessible. 
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